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Monte Carlo methods. Figure ^ shows the geometry 
used in the "2.45 MeV w/shielding" detector response 
simulation. The cavitation fluid, modeled as described in 
footnote 13 of I'l, consists of carbon, deuterium, chlorine, 
oxygen, nitrogen, and uranium. The simulation, using 
Geant4 includes all relevant neutron interactions, 
particularly elastic scattering, thermal capture, (n, n'7) 
scattering, and neutron-induced fission. 

To calculate the response function, neutrons of energy 
2.45 MeV, emitted isotropically from the center of the 
flask, scatter through the materials. When a neutron 
elastically scatters protons in the liquid scintillator, the 
recoil energies are converted to equivalent electron en- 
ergies 0, summed, and then smeared according to the 
detector's resolution function, eventually obtaining the 
response function g{Ei,e)- 

In the same manner, I calculate the other response 
functions with the following changes: the "2.45 MeV" 
simulation does not include the paraffin shield, and the 
radioisotope 0,0 simulations "Cf-252" and "PuBc" as- 
sume there are no intervening scattering materials be- 
tween the sources and the detector. 

Statistical methods. Following the notation of [Qj , the 
raw data from Fig. 9(b) of ^ are 

rifc = "cavitation off" counts in channel k 
n'l. = "cavitation on" counts in channel k. 



Each run is 300 s in duration, and Fig. 4 of shows the 
background-subtracted signal n'j. — rik- 

The Uk background data are modeled by a sum of two 
exponentials, and the n'j, data are modeled by the same 
background function plus the scaled response function. 



yk = Ai exp(-fc/yl2) 
y'k = Vk + Asgk. 



A3 exp(-/c/^4) 



The binned response function gk, is found by averaging 
g{Eee) over the energy range of channel k. 
Then, the Poisson likelihood chi-square 6j is 
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where 
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k = Uk - nk + nk\n{nk/yk) 
'^'k = y'k-n'k + n'k\TL{n'k/y'k). 



Note that, under proper conditions '6'|, asymptoti- 
cally approaches a distribution. Moreover, better fits 
give lower values of Xa p- Minimization of Xa p deter- 
mines the five fit parameters At. See Fig.|3[a) for the fit 
using the "Cf-252" response function. 



To determine the distribution /(xa p) f^'' ^ given fit, I 
sample from many synthetic data sets, each chosen, for 
k — W, . . . , 249, from Poisson distributions of mean value 
yk and y'f, . In the Comment, I report the goodness-of-fit 
as a Z-value, defined by 



obs. 



which expresses the observed value of Xa p i'^ terms of 
the equivalent number of standard deviations from the 
mean of a normal distribution. As shown in Fig. I^^b), 
the observed value of Xa p foi' the "Cf-252" fit is within 
one equivalent standard deviation and is therefore sta- 
tistically consistent. The other three fits are outside five 
equivalent standard deviations, and are therefore statis- 
tically inconsistent. 
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FIG. 2; Fit statistics, (a) Fit using the simulated ^^^Cf response function for g^. The histograms are and nj,, and the 
smooth blue lines are theoretical curves and y'f.. The minimized value of Xa,p is 432. For clarity, the two graphs are offset 
by ten counts, (b) Numerically sampled distributions of Xa,p for the four hypotheses. Observed values of Xa,p ^^re shown. 



